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Very recent experiments have studied for the first time 
collective excitations of an ultracold 6 Li gas [1,2], cov- 
ering in particular the BEC-BCS crossover domain [2]. 
We wish to point out that the results for the axial mode, 
through hydrodynamics, give direct access to the (3D) 
equation of state of the strongly interacting gas, mostly 
near the unitarity limit. On the other hand the surpris- 
ing results found for the radial mode are actually not 
necessarily in contradiction with the expectations from 
superfluid hydrodynamics. 

Indeed the radial mode frequencies tt r are near [1] 
2tt.2400 Hz and [2] 2tt.1200 Hz. The frequency uj f = 
Ep/Ti corresponding to the Fermi energy is [2] 1.6 10 5 
Hz. If on the BCS side we estimate the gap by the 
BCS relation A = 1.75T C and take a nearly optimal value 
T c = 0.2Ep, we have M7 r /A ~ 0.13 at the trap center, 
this ratio increasing when we go away from the center. 
This is large enough to make questionable hydrodynam- 
ics, which assumes M7 r /A <§; 1 in order to be accurate, 
and may explain the 10% discrepancy between theory 
and experiment at unitarity. Similarly one can estimate 
that the ratio l;/l± between the Cooper pair size and the 
transverse size of the trap is at best given by this same 
figure ~ 0.1, which makes doubtful the accuracy of the lo- 
cal density approximation (a necessary ingredient in the 
hydrodynamic result). The situation in Rcf. [1] is even 
worse. This is also consistent with a natural superfluid 
interpretation of the strong attenuation [2] at 910 G as 
a pair-breaking peak corresponding to hfl r = 2A(T, £?), 
occuring because T c and A decrease with increasing field 
B. Clearly in this case KCl r / A is no longer small. Taking 
for example A(T, B) ~ T c (the gap depends on T/T c , 
which is not known) gives T c /Ep ~ 0.02, coherent with 
the estimated temperature T in this experiment. 

On the other hand the axial mode is a very good case 
for making use of the hydrodynamic limit. Indeed its fre- 
quency is very low hfl a /A ~ 5.10~ 3 , which is the appro- 
priate range for this approximation. Moreover the exper- 
imental temperature is certainly quite low, as confirmed 
by the very low damping found in most of the magnetic 
field range. This makes it possible to neglect the effect of 
dissipation on the frequency of the modes. Finally for the 
very elongated traps used in experiments, one has to deal 
with a simple effective one-dimensional problem. In this 
case, when the chemical potential /i(n) of the (3D) gas is 
a power law of the total particle density n 1 ^, the axial 
mode frequency is given by an exact analytical result [3,4] 



u> 2 /uj 2 = 2+l/(p+l). In the present situation this case is 
found in the BEC limit (small positive scattering length 
a), with p — 1, leading to uj 2 /uj 2 = 5/2, in the BCS limit 
(small negative scattering length a), with p = 3/2 and 
uj 2 /uj 2 = 12/5. Moreover quite remarkably [5] this same 
value p = 3/2 applies also in the unitarity limit, where 
a is very large. These values are in fair agreement with 
experiment [2] for the BEC and the unitarity cases. 

We have shown recently [6] how it is possible to invert 
such experimental data to obtain the equation of state of 
the gas. This requires only that fJ.(n) is known in some 
limiting case, from which one then go away by an iterative 
procedure, making use of the experimental knowledge of 
the mode frequency as a function of density. The basic 
ingredient of this method has been shown recently [7] to 
have an accuracy of order 10~ 3 . In the present case we 
have in principle the choice between the three limiting 
cases mentionned above, the most convenient one being 
the unitarity limit. However there are not enough data 
to carry out the above program with a sensible precision. 
Nevertheless we can analyze the region in the vicinity of 
the unitarity limit, where a fairly linear behaviour is ob- 
tained experimentally. In this case we have the simpler 
problem of performing a perturbative calculation [71 . 
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FIG. 1. Reduced axial mode frequency as a function of the 
inverse scattering length a -1 for the model in the text. Heavy 
line: linear approximation near unitarity. Red squares: ex- 
perimental data of Ref. 2. 

In the vicinity of the unitarity limit we linearize the 
general form /u(n) = h 2 'k 2 F } '(1 / 'fcpa) /2m of the chemical 
potential by f(y) = £ — Sy, where kp = 37r 2 n and £ and 
S are constants. Using the method of Ref. [7], we find 
[8] that, in this regime, the shift of the mode frequency 
is given analytically by S(uj 2 /uj 2 ) = uj 2 /uj 2 — 12/5 = 



1 



(256/875Tr)(S/£)(l/k Fmax a) where k Fmax is the (3D) 
Fermi wavevector at the center of the trap. Taking [2] 
E F = 1.2/j.K and 5(1/ a) = 5.7 itfm^G- 1 in the vicin- 
ity of the resonance, this gives 5(lu 2 /u>^) ~ 10 _3 S/£G _1 . 
Comparing with the experimental result ~ 1.1 10~ 3 G _1 
we obtain S/£ = 1.1. If we take for £ = 1 + [3 a value 
[5] £ = 0.45 which is most likely both experimentally 
[9] and theoretically [10], we find the experimental result 
S~0.5. 

Finally a quite simple model in reasonable agree- 
ment with known constraints is f(y) = 1/2 — 
(1/tt) arctan(7ry/2). Its linearization for y ~ leads 
to take £ = S = 0.5. It gives the proper limit for the 
weakly interacting Fermi gas and in the BEC limit it 
yields a m — 6a /tt for the molecular scattering length, 
which is not too different from the result a m — 0.6a of 
Petrov et al [11]. Although it is very simple this model 
gives already a quite good agreement [12] with experi- 
ment [2] in the resonance region as can be seen in Fig.l. 
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